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We study the scattering of a matter-wave from an interacting system of bosons in an optical lattice, focusing 
on the strong-interaction regime. Analytical expressions for the many-body scattering cross section are derived 
from a strong-coupling expansion and a site-decoupling mean-field approximation, and compared to numerically 
obtained exact results. In the thermodynamic limit, we find a non-vanishing inelastic cross section throughout 
the Mott insulating regime, which decays quadratically as a function of the boson-boson interaction. 
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I. INTRODUCTION 


In a many-particle system, complexity owed to its multi¬ 
partite nature arises from two fundamental sources: the in- 
distinguishability of the particles constituting the system and 
their mutual interactions. In today’s ultracold atoms exper¬ 
iments, ensembles of interacting identical particles can be 
trapped in optical potentials |[Il,|2l and controlled with extreme 
precision, offering the opportunity to study complex many- 
body dynamics beyond mere single-particle or mean-held de¬ 
scriptions. Simulations of solid-state and strongly-correlated 
physics have come into reach Arguably most prominent 

was the observation of the superhuid to Mott insulator phase 
transition of bosons in an optical lattice ||3l, but also Fermi- 
Hubbard models, in and out of equilibrium ifioljisl] . 

have been realised. Currently, the effects of disorder lll4l - [l^ 
and its interplay with interactions iflol - E^ are objects of inten¬ 
sive research, as are implementations of laser induc ed g auge- 
helds for the simulation of quantum magnetism 

One challenge that emerges in these systems concerns the 
extraction of the relevant information. While standard time- 
of-hight (ToF) measurements allow to infer the expansion 
dynamics of the system Q, more elaborate schemes like 
Bragg-spectroscopv permit to extract, for instance, spectral 
information Whereas ToF-experiments require the 

destruction of the sample in each run of the measurement, 
direct imaging of the s}^tem can be achieved in-situ by the 
scattering of light ll4^5in. w ith current experiments reaching 
single-site resolution 11521 - 1^ . 

Inspired by neutron scattering j^lbll] . the scattering of an 
atomic matter-wave has recently been piU forward as a com¬ 
plementary in-situ imaging technique 11621 - 1^ . Additionally, 
matter-waves in combination with cold-atom crystals appear 
as convenient simulators of low-dimensional solid-state sys¬ 
tems ibSl - l^ . Experimentally, a BEC has been used as probe 
of a cloud of bosonic atoms in an optical lattice with the 
elastic scattering signal reflecting the crystal structure of the 
strongly-interacting system, while, in the weakly-interacting 
regime, inelastic scattering induces single-particle band exci¬ 
tations. It has been shown in earlier works that the inelastic 
scattering signal of the matter-wave probe decays as a function 
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of the strength of the interaction between the target bosons 
1^ P^ . In the regime where this interaction strength 
is too weak to deplete the condensate significantly, we have 
successfully employed Bogoliubov’s theory to describe ana¬ 
lytically the scattering cross section ll^P^ . There, we found 
that the decay of the inelastic signal always emerges linearly 
and is independent of the system size, number of particles, 
and condensate depletion. 

In this work, we develop an analytical description of the in¬ 
elastic scattering cross section for strong interactions. In par¬ 
ticular, we address the question of whether the inelastic cross 
section vanishes at the superfluid to Mott insulator phase tran¬ 
sition. In Sec. m the stage is set by introducing the cross sec¬ 
tion of a many-body target. In Sec. |ml analytical expressions 
for the inelastic cross section for strong interactions are de¬ 
rived from a strong-coupling expansion iSec. lIII All , and from 
a mean-held approximation [Sec. IIIIDII . In Sec. IIVI we per¬ 
form a thorough analysis of our findings, where we compare 
the analytical predictions obtained in Sec.|nI]to exact numer¬ 
ical simulations, before we conclude in Sec.|V] A few rather 
technical derivations are left to the appendices. 


II. MANY-BODY SCATTERING CROSS SECTION 
A. The Bose-Hubbard model 


The target that we consider is a collection of interacting 
bosons in the lowest single-particle energy band of a one¬ 
dimensional optical lattice with L sites, described by the 
Bose-Hubbard model lEH: 



1) — gLTlj 


(iJ) 


where J > 0 is the strength of the nearest-neighbor hop¬ 
ping, given in the latter Hamiltonian by the sum over all pairs 
{i,j) of nearest-neighboring sites. The onsite boson-boson in¬ 
teraction strength U is derived from a Eermi-pseudopotential 
if^ . and we restrict our considerations to repulsive interac¬ 
tions with C/ > 0. In the above grand-canonical description, 
the chemical potential jj, Axes the average number N of parti¬ 
cles. The operators ay and dj are respectively the annihilation 
and creation operators of a boson on site j, obeying bosonic 
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commutation relations, 

[ai,al] = Sij, [aj,aj] = [aj,a]] = 0, (2) 

and fij = aJjQj counts the number of particles on site j. 
The tunneling strength J and the interaction energy U can 
be controlled by tuning the optical lattice and by means of 
magnetic Feshbach resonances IiIItI. Throughout this 
work, all energies are measured in units of the recoil energy 
Er = h'^k‘l/2M, where M is the mass of one boson and k]^ is 
the wavenumber of the laser providing the lattice. The natural 
length scale is given by the lattice constant d = tt/Zcl- 

Depending on the value of U/ J, the Bose-Hubbard model 
exhibits two distinct thermodynamic phases s. In the non¬ 
interacting limit where U/J 0, the ground state of JTbh 
is a compressible superfluid (SF) state, in which all bosons 
condense into the same delocalized single-particle Bloch state 
of the lattice. In the opposite limit U/J oo, the bosons 
localize on the wells of the lattice, and for an integer filling 
factor 


N 

(3) 

the ground state of K^h is given by a Fock state with exactly ly 
bosons on each site, called a Mott insulator (MI). In the ther¬ 
modynamic limit, the system undergoes a quantum phase tran¬ 
sition from the SF to the MI regime for a finite value ofU/J. 
The SF phase is characterized by a non-vanishing superfluid 
fraction, compressibility and on-site number fluctuations, and 
a vanishing excitation gap. In contrast, at the transition point 
to the MI phase an energy gap opens, occupation number fluc¬ 
tuations are suppressed due to the competition between kinetic 
and interaction energies, the system becomes incompressible 
and the superfluid fraction vanishes 


B. Scattering from a Bose-Hubbard target 

We use the formalism presented in Refs, pT^ . which 
for the benefit of the reader we summarize in the following. 
We study the scattering of an atom (probe) of mass m from 
a Bose-Hubbard system (target), and assume that the optical 
lattice is transparent to the probe. Such a species-selective 
optical lattice —Phased on a species-selective dipole potential 
117^ — has already been exp erimentally realized in a mixture 
of ®^Rb and atoms ||^|^. The incoming energy of the 
probe is considered to be such that no interband excitations 
of the system can occur, and s-wave scattering dominates the 
interaction with each target atom. The interaction can there¬ 
fore be described by the s-wave scattering length as through 
a pseudopotential: 



FIG. 1. (Color online) Scattering setup: A particle of mass m ini¬ 
tially in a plane-wave state with momentum fein is scattered into the 
angle 6 from a target of atoms (all of which have mass M) submerged 
in a one-dimensional optical lattice with depth Vo and lattice constant 
d = Tr/fer., where /cl is the laser wavenumber. The asymptotic final 
state of the probe has momentum k. Figure taken from Ref. j^. 


probe’s asymptotic incoming state is a plane wave |fein) with 
wavevector fcin, and incoming energy Ein = h^kf^/2m. In 
Born approximation, the far-held scattering cross section is 
given by 


1 da _ A — Eq 

Ellm~ e/. 


dr {(j)n\n{r) \(j)o) 


( 5 ) 


where {\4>n)} are the eigenstates of Kbh with corresponding 
eigenenergies {En}, and k is the transferred momentum: 

K = fein - fe, (6) 


with fe the asymptotic outgoing momentum of the probe. We 
choose the lattice axis along the Ma;-direction, thus Tj is given 
by Tj = XjUx = jdux- For normal incidence of the probe 
with respect to the lattice axis, the scattering is invariant under 
rotations around the x-axis (neglecting interference terms be¬ 
tween the incoming and scattered waves), and we can restrict 
our considerations to the x-y plane, as shown in Fig. □ The 
transferred momentum along the direction of the lattice 


K = K- Ux, 


( 7 ) 


obeys 


where 


K = Kel 



(8) 


Keic? = —7rsin0 


m Ei„ 

M Ex 


( 9 ) 


V{r) = ‘^^^as^S{r - (4) 

™ p=i 

where r and {'r(^A/3=i....,Af give the positions of the probe 
and the target atoms, respectively. The target is assumed to 
be initialized in its ground state |())o), with energy Eq, and the 


is the X component of the transferred momentum for elastic 
scattering. 

The bosonic field operators 'ijj{r) and annihilating 

and creating a particle at position r, respectively, define the 
density operator 


fiir) = 


( 10 ) 
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and they can be expanded in the single-particle Wannier basis; 

L 

='^ajw{r-rj), ( 11 ) 

i=i 

where w{r — Vj) is the Wannier function describing a particle 
localized around the position of the j-th lattice site 11^ 1^ . 
From expansion (fTTl i the density operator is split into diago¬ 
nal and off-diagonal contributions in the Wannier basis, the 
latter being proportional to the overlap of Wannier functions 
centered at different sites. We consider a deep lattice of depth 
Vh = 15Er, for which this overlap can safely be neglected 
11^ . In this diagonal approximation, the inelastic cross sec¬ 
tion (i.e. when the target is left in an excited state) reads 


1 da 




En — En 


Eir 


i=i 


\WiK)\ , (12) 


where W (k) is the form factor of a unit cell of the lattice: 


W{k) 




( 13 ) 


The inelastic cross section bears a clear signature of the na¬ 
ture of the ground state. In the SF regime, fluctuations of the 
on-site occupation numbers give rise to a non-vanishing in¬ 
elastic cross section, which can be described a naly tically both 
at [/ = 0 (from the exact solution of Kbh) II^ . as well as 
in the regime of weak-depletion of the macroscopically occu¬ 
pied single-particle ground state 1^ . using Bogoliubov’s the¬ 
ory. The inelastic cross section decays as a function of U/J, 
and in the MI limit (Uj J oo), it vanishes exactly in the di¬ 

agonal approximation (see Table I in Ref. as the ground 
state becomes an eigenstate of the on-site number operator fip 
In this work, we want to reveal how the inelastic cross sec¬ 
tion depends on the interaction in the regime of strong, but 
finite interactions. To this aim, we employ a strong-coupling 
expansion (SCE) of Hamiltonian ([T]), in which a perturbative 
series of its eigenstates and energies is constructed by regard¬ 
ing the tunneling operator as a perturbation to the pure inter¬ 
action Hamiltonian Moreover, we also show how a 

site-decoupling mean-field approach fails to describe the be¬ 
havior of the cross section. 


III. DERIVATION OF THE INELASTIC CROSS SECTION 
IN THE STRONGLY-INTERACTING REGIME 


A. Strong-Coupling expansion 


In order to set up a perturbative expansion, the grand- 
canonical Bose-Hubbard Hamiltonian [Eq. ([T])] is split into 
the unperturbed on-site part. 


ifo = E 

j=i 





and a perturbation given by 

L 

V = J2 , ( 15 ) 

i=i 

with all energies in units of U, indicated by the tilde notation, 
i.e. fl = p,/U. Eor small values of the dimensionless tunneling 
strength J = J/U, the eigenstates and eigenenergies of the 
full Hamiltonian 


Kbh = Ko-JV ( 16 ) 


can be obtained as a power series in J. Notice that iCsH, 
Kq and V, all commute with the total number operator N = 
J2f=i therefore all energy eigenstates can be chosen 

to have a well-defined number of bosons. 

The eigenstates of Kq can be taken as eigenstates of the 
on-site number operators fij. 


|n) = |ni,...,ni), 
fij |n) = Uj |n), Uj e No, 


with unperturbed energy 8^'- 

Kq |n) = |n), 


L r 


£n=E 

i=i 


- 1 ) - jlrij 


( 17 ) 

( 18 ) 


( 19 ) 

( 20 ) 


Eor integer filling factor v [Eq. Q], the ground state of Kq, 

|0mi), (21) 


is non-degenerate. It corresponds to z/ particles localized on 
each site, and has unperturbed energy 


— L 




( 22 ) 


We write the perturbation series of the ground state |(^o) of 
iTfin as 


|<Ao) = ^0 |<^o), (23) 

\4>Q) = \K^-h4>^Q'') + 0{P), (24) 

Eo = - Je''^^ + + 0(P). (25) 

To first order in the tunneling, the correction to this state 
is readily obtained from non-degenerate perturbation theory 
(PT): 


\9q )= 2 ^ p ( 26 ) 

Hence, the first correction is a superposition of states with one 
additional particle-hole (PH) pair, i.e. of states of the form 
|n) = liy,... , 1 / ± ^ 1,... , 1 /), with excitation energy 

F — F ■ — ^ 

<-'mi — -L* 


( 14 ) 


( 27 ) 














4 


In the ground state, the first-order correction to the energy 
= (</. mi I V \ 4>mi) vanishes, since V only connects states 
that differ by one PH pair. More generally, all odd-order cor¬ 
rections to the ground-state energy can be shown to vanish. 
To third order, the perturbed energy of the ground state is ob¬ 
tained straightforwardly HH] ; 

Eo=£mi-2Liy(iy+l)J^ + 0(J^), (28) 

Symbolical expansions for ground-state properties can be ob¬ 
tained on a computer within non-degenerate perturbation the¬ 
ory to very high orders ll^ [ssl ■ 

The inelastic cross section [Eq. (fl^ l is governed by the 
matrix elements 


{(l)n\nj\(j)o), (29) 

which are non-zero only if the states |(/)„) and \(j)Q) have the 
same number of particles. Therefore, the SCE of ( l29l l requires 
also the perturbative series of the fixed-density excitations of 
Kq, which can be characterized by the number of PH pairs. 
Given that the tunneling operator changes the number of PH 
pairs by one in every subsequent term of the perturbation se¬ 
ries of a state, the desired order in J for the matrix elements 
dictates which manifolds of excitations need to be taken into 
account: Eor the expansions of ( |29] | to order r, the perturba¬ 
tive series to that order of all iTo-excitations including up to 
r PH pairs, and the ground state, must be considered. To first 
order, only one-PH states of the form 


|s, 5 -f /) 


alas+i 
\/v{v -b 1) 


I'/’mi) ) 


s = 1,... ,L, 

Z = 1,...,L- 1, 


(30) 


contribute, where s and s + l denote the position of the particle 
and the hole, respectively, and periodic boundary conditions 
are used JH. The restriction on the distance I enforces that 
particle and hole be on different sites. The states (l30l l span 
an L{L — l)-fold degenerate subspace V with unperturbed 
excitation energy 


^ph ^mi — 1- 


( 31 ) 


Thus, we need to employ degenerate PT, and find states with 
a well-defined limit as J —0. The latter are obtained by 
diagonalizing the perturbation V in the degenerate subspace 
22. Translational invariance of K^h allows to separate 22 into 
L manifolds characterized by different well-defined center-of- 
mass (CoM) quasimomenta q iF^[83.18^ : 




q = 2'Kj/Ld, 
j=0,...,L-l. 


(32) 

In this basis, the CoM degree of freedom decouples and the 
perturbation becomes block-diagonal 1^ : 

{q', l'\V\q, 1) = + T,Sv,i+i] , (33) 

with the effective (complex) tunneling parameter 

Tg = {iy + -b (34) 


The matrix (l33l l describes the problem of a particle on an open 
chain with hard-wall boundary conditions at sites L and a fic¬ 
titious site 0 (since I has to be different from both, 0 and L). 
Diagonalization is achieved by the states 


L-l 


9, >^) = y 7 E sin(>rd k: 0 > 


= TTi'/Ld, 

1=1 f = l,...,L-l, 

(35) 

where we note that can be interpreted as the quasimomen¬ 
tum connected to the PH relative coordinate. The phase (q 
is the argument of the complex tunneling parameter Tq, and 
given by 


tanG = 


{i' -b 1) sin/gd) 


v + {ly + 1) cos{qd) 

The eigenvalue of the perturbation (restricted to 22) is 


(36) 


= 2 cos{>cd)\/l + -b 1) cos2(gd/2), (37) 


which gives the first order energy correction of the eigenstates 
\4>q,>e) of Kbh whose (J = 0)-limit are the one-PH states 


= |g, - J |</.W ) + 0(j2), (38) 

Eq,^ = £ph - (39) 

Note that to first order the energy is still degenerate in the 
center-of-mass quasimomentum q. This degeneracy, however, 
plays no role for the perturbative expansion, since, as men¬ 
tioned above, the perturbation is translationally invariant and 
thus never couples states with different values of q. One can 
then proceed as if the degeneracy were fully lifted i^ . The 
second order correction to the energy is given in AppendixlAl 
The first-order correction to the state can be written as 


kgli) = (kmil^k; l</>mi) “f kc) + \vv) ■ (40) 

Here, the first term corresponds to the ground-state contri¬ 
bution, \vc) to a state orthogonal to both, |(/>mi) and 22, and 
\vT>) is the contribution of all states from within 22 with 
^ X. The explicit expressions to obtain \vc) and jui)) 
are given in AppendixlAl However, for a first-order expansion 
of the matrix element (l29l l. only the ground-state contribution 
in Eq. ( 00 ]) is necessary, since \vc) and |t;x>) couple to ^o) 
[Eq. (l24l) 1 only to higher orders in J. 

We would like to emphasize that SCE provides a rather ac¬ 
curate description of the phase diagram of the Bose-Hubbard 
model, which can be obtained from the energy gap between 
the ground state and the lowest non-number-conserving exci¬ 
tations of the system (i.e. adding/removing one particle) lHH- 
ImI] . The resulting equation for the critical point of the SF-MI 
phase transition at fixed density is given in AppendixlAl 


B. Inelastic cross section from SCE 

The leading contribution in J of the inelastic cross section 
[Eq. ([T2l) 1 is obtained by considering all matrix elements (l29b 
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with a non-vanishing first order term. As explained above, 
only the manifold of eigenstates \(t>q, 3 <) needs to be taken into 
account. From the expansions (l24l) and d^ . the matrix ele¬ 
ment of the onsite-density operator is found to be 

{4>q,^\nMo) = J {q, x\[nj,V]\4':ni) +0(J^) 

= + 0(P), 

( 41 ) 


with M{q, k) given by 
M{q, k) = — 2i sin(>r(i) sin(gd/2) 


(42) 


Inserting Eqs. dMll, <l39ll and (gB into Eq. ([I2ll, we obtain an 
expression for the inelastic scattering cross section, valid for 
large values of the interaction strength; 


1 da 

SCE 

■ J' 


Nal dn 

inel 

_u_ 

/-V Bin 

q,>c ’ 


X 


M{q, >c)T,{k 


2 


q)W{K) 


(43) 


where the excitation energy gap has been evaluated to first or¬ 
der in J only, and we normalized to the number of bosons N, 
the scale of the inelastic cross section in the non-interacting 
limit 1^ . The function — q) = rep¬ 

resents the interference of elementary waves emanating from 
each scattering centre. 


|S(«-g)p 


sin^ ((k — q)dL/2) 
sin^ ((k — q)d/2'j 


(44) 


Notice that, although not indicated explicitly, the transferred 
momentum k [Eq. ®] depends on the transferred energy, and 
thus on the quasimomenta q and via 


K = Kell/l - 


U{1 - JE^qll) 
E\„ 


(45) 


Equation (l43l l shows that, for large values of the interaction, 
the inelastic cross section decays quadratically with U/ J, in 
agreement with the behavior of the dynamic structure factor 
of a Bose-Hubbard system j^ . The dynamic structure factor 
is intimately linked to the scattering cross section, being the 
response of the system to a density perturbation 


C. Site-decoupling mean-field description 


where in one dimension the mean-field coupling strength 
reads 


Xj — J^(aj+i) + (dj-i)'^, (47) 

and the expectation values are taken with respect to the 
ground state of the system. In a translationally invariant lat¬ 
tice, and due to the invariance of ATbh under a global U(l)- 
transformation, the above quantity can be chosen to be real 
and constant for all lattice sites. Then, the ME coupling 
strength A = Xj is fixed by the self-consistency relation 

X = 2J{dj), (48) 

since the expectation value itself depends on A. The resulting 
effective Hamiltonian is a sum of single-site operators: 


ffeff = = 


U) 




j=i 


i=i 


(49) 


where 

^>-0^ = hj, (50) 

y(o) =-(^a] + dj), (51) 

and A = A/C/. Numerically, it is found iteratively that A is 
the mean-field order parameter of the SE-MI phase transition 
iItM Isol - loil] : it is zero in the MI phase, and takes on non¬ 
zero values continuously as the phase boundary into the SE 
regime is crossed. Therefore, in the vicinity of the critical 
point, A is small, which motivates a perturbative expansion of 
the eigenstates and eigenenergies of iFeff in A. The fact that 
Hgff is a sum of single-site operators allows us to reduce the 
calculation to the single site Hamiltonian h = . 

The single-site eigenstates \v) of /ip are characterized by 
the integer number of bosons v, and have unperturbed ener¬ 
gies 


- 1) - (52) 

in units of the interaction energy U. In order for |z/) to be the 
state of lowest energy, the chemical potential has to satisfy 

V -1 < fi <v, z/eN. (53) 

The perturbative expansion of a generic eigenstate |^i,) of h 
reads 


k.) = + A + A^ + OCX% (54) 


An effective Hamiltonian results from ([T]) when the tunnel¬ 
ing (the coupling of a given site to its neig hbors) is treated in 
terms of a mean-field (ME) coupling iItM l89l - l^ . i.e. when 
the full hopping term in ([Til is replaced by 

L 

- JX] aloj — >-Yl 

(ij) 3 = 1 


where the corrections follow from standard non-degenerate 
perturbation theory: 

k - 1) - + 1): (55) 

\v-2)+ \v + 2) (56) 

+ ( 57 ) 
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in terms of the coefficients 


regime: 


(fc) ^ + [1 -sgnA:]/2 

e{v) — £{v + k) 


k^O. 


(58) 


Similarly, the perturbative expansion of the energy ii, of |^i,) 
reads 


1 da 

Na1 dfl 


MF 

inel 


A2 


1/ 


E 



E\„ 




( 66 ) 


with the transferred momenta 


e, = e{v) + + vVTTcW) P + ^(A^). (59) 

Using the second-order expansion of the value of the 
MF coupling strength A in the SFphase is estimated from the 
self-consistency condition ( l48l l 11^ : 

A + Vi'd- (61) 

B =\/u + 2 -f y/v -I 

D. Inelastic cross section from the MF description 


= «ei-\/EcEEE, (67) 

and it must be emphasized that A^ = 0 for J ^ , and A^ 

is given by Eq. (l60l) for J > where 

jff =1/2 + V- ^v{v + l) (68) 

corresponds to the MF estimate of the FD critical point. 
Hence, according to the MF description, the inelastic cross 
section should decrease linearly with [// J in the vicinity of 
the transition when approached from the SF phase, vanish at 
the critical point and remain zero in the whole Mott phase. As 
one goes deeper into the SF phase (and thus the value of A 
increases), a higher order expansion of the single-site states 
would be necessary to estimate correctly the dependence of A 
on [7/ J, which will include non-linear terms. 


Due to the separability of lEq.(l49ll1. the calculation of 
the matrix elements (|2^ in the inelastic cross section reduces 
to the single-site level. We therefore need to evaluate the ma¬ 
trix element for the ground and excited states |^i,) 

and respectively, of the single-site Hamiltonian h. From 
Eq. ( |55] | one can see that, to first order in A, only excited states 
with v' = V ± \ must be considered: 

(^,+i|n|e.) = -AcW+0(P), 

(^.-i|n|e.) = Ac[,-^)+0(P). (63) 

The contributing many-site excited states are of the form 
\^u±i)j ICiz)*, for j = 1,..., L. Therefore, to leading 
order in A, the inelastic cross section [Eq. (IT^ I is given by L 
times the single-site contribution if^ . To first order in A, the 
energy gap to the contributing excited states is given by [see 
Eq. ®] 

Af =£{v±l)-£{v). (64) 

Finally, to derive the MF cross section, the value of the 
chemical potential fi needs to be specified. We aim at a de¬ 
scription of the scattering cross section in the vicinity of the 
fixed-density (FD) SF-Ml transition. The corresponding value 
of the chemical potential can be obtained from the MF phase 
diagram of the Bose-Hubbard model iFt^ (for details see, e.g.. 
Ref. izi): 

-f 1) - 1. (65) 

Putting everything together yields the MF inelastic scattering 
cross section to leading order in A, in the strong interaction 


IV. ANALYSIS OF THE INELASTIC CROSS SECTION 

In order to gauge the quality of the approximate analytical 
expressions (14^ and (i66l l for the inelastic cross section, we 
calculate numerically the exact cross section, Eq. (IT^ . via ex¬ 
act numerical diagonalization of Kbh in a system with fixed 
integer density. The Wannier functions are approximated by 
Gaussians, corresponding to the ground-state wave function 
of a harmonic approximation of each optical potential well, 
leading to a Gaussian form factor [Eq. mi 

WiK) = e (69) 

This is a good and common description for deep enough op¬ 
tical lattices, and valid for our choice of the lattice depth 
Vq = 15Er, which corresponds to the tunneling strength 
J = 6.5 X 10~^Er, and gives rise to a gap to the second 
band of 6.28Er 1^ . 

A. ME expression 

Figure|2]shows the mean-field prediction [Eq. (i66l) 1 in com¬ 
parison to exact numerical results for different system sizes 
at unit filling. Clearly, the mean-field approximation fails to 
describe the behavior of the inelastic cross section. As dis¬ 
cussed in Sec. lIIIDl the MF cross section exactly vanishes in 
the whole Mott phase up to the estimated critical point 
[Eq. (i68l l1. Comparison to numerical results, however, show 
that—at least for finite system size— the inelastic cross sec¬ 
tion is non-zero in this region. For the system sizes acces¬ 
sible in our numerical simulations, no trend towards the L- 
independent MF result (i66l l can be observed. 
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FIG. 2. (Color online) Inelastic scattering cross section in the strong- 
interaction regime for unit filling factor i/ = 1. The symbols show 
exact numerical results for different system sizes, while the red solid 
line corresponds to the analytical formula J66b . obtained from MFT 
(only in the region where A ^ 0.1). The latter predicts the vanishing 
of the inelastic cross section at the critical interaction strength — 
11.66J, indicated by a vertical dotted line. Relevant parameters are 
Ei„ — 2Er, m = M, and 0 — 0.99. 


A quantitative disagreement between the simulations and 
the MF result should be expected, given the known poor per¬ 
formance of MF in one dimension. More elaborate mean- 
field like approaches (e.g., a dynamical MF approach lH^l) 
might provide a better description of the inelastic cross sec¬ 
tion. However, the discrepancy found here clearly demon¬ 
strates that the site-decoupling MF formalism does not yield 
a qualitatively correct result in one dimension, if employed in 
the above straightforward manner. 



FIG. 3. (Color online) Upper panel: Inelastic scattering cross section 
for a system with L = 4 sites and different bosonic densities i/. For 
large enough interaction U/J, the SCE formula l l43b (lines) correctly 
describes the exact numerical data (symbols). Vertical dotted lines 
indicate the critical value of the SF-MI phase transition for each 
V {v = 1, 2, 5 from left to right), as predicted from third order SCE 
(see Appendix 13. Lower panel: Relative difference Aics [Eq. l l70l ll 
between the analytical approximation and the exact numerics, versus 
the interaction strength in units of the respective critical interaction, 
for different system sizes and densities. The horizontal dotted line 
indicates a deviation of 10%. The non-monotonicity of Aics visible 
for 1/ = 4 is a finite-size effect due to a crossing of the analytical 
approximation and the exact numerical data. The inset shows data 
for L ^ 6 on a semi-logarithmic scale. Relevant parameters for both 
panels are E^^ = 2Er, m = M, Q = 0.99. 


B. SCE expression 

In contrast to the simplistic site-separability of the eigen¬ 
states of the system enforced by the MF approach for all val¬ 
ues of U/J, the SCE of Sec. IIII Al describes the eigenstates 
of iCBH for J 0 in terms of superpositions of Fock states 
with different occupation distributions, which are in general 
non-separable in the site basis. 

The SCE treatment leads to the L-dependent expression 
(l43l l for the inelastic cross section, which is predicted to de¬ 
cay quadratically for large U/J, and hence to be finite —at 
least for finite system size— inside the Mott phase. This re¬ 
sult is confirmed in the upper panel of Eig. [3 where Eq. (l43l l 
is compared to the numerically calculated exact cross section 
for a system size of L = 4 sites and different filling factors 
V. It is however apparent that the higher the density, the larger 
the interaction has to be in order for the analytical approxima¬ 
tion to be valid. This is correlated with the fact that for higher 
densities the critical point of the SE-MI transition is shifted to¬ 
wards stronger interactions S (in the upper panel of Eig. [3 
the corresponding positions of the critical points ji''\ esti¬ 
mated from SCE as described in Appendix lAl are marked by 
vertical dashed lines). This behavior seems to indicate that 
the parameter controlling the validity of the perturbative ex¬ 


pansion is not J but rather the ratio J/ J^'^. This conjecture 
is further supported by the lower panel of Eig. [3 which shows 
the relative difference 


da 

exact 

da 

SCE 

dCl 

inel 

dO. 

inel 


da 

1 exact 



d^l 

inel 



between the exact inelastic cross section [Eq. ([T2li l and the 
SCE expression [Eq. (143) 1. as a function of U in units of 
the critical interaction Uc''\ for a fixed scattering angle. In 
terms of the renormalized interaction strength, the relative 
deviations Ajcs for constant L and different v nearly coin¬ 
cide, and they exhibit the same behavior when L is increased. 
Eor U > Uc'^'^ the Mott phase is approached and the devia¬ 
tions become smaller with increasing interaction, whereas for 
U < Uc'^'^ the deviations become increasingly pronounced. 
The non-monotonicity of Aics visible for L = 4 is due to 
a crossing of the analytical approximation and the exact nu¬ 
merical data; a finite-size effect which disappears for L > 5. 
The inset of the lower panel of Eig. [3 reveals that, around the 
critical point, Aics depends exponentially on the renormal¬ 
ized interaction strength, with an exponent that changes as the 
transition is crossed and seemingly converges to a system size 
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FIG. 4. (Color online) Inelastic scattering cross section in the strong- 
interaction regime for unit filling factor v = 1. The symbols show 
exact numerical results for different system sizes, while the thick 
solid line corresponds to the L ^ oo limit [Eq. l l72b l of expres¬ 
sion (Ell, obtained from SCE. The vertical dotted line indicates 
the critical point of the SE-MI phase transition obtained from third- 
order SCE (see Appendix lAt. Relevant parameters are Em = 2Er, 
m = M and 6 = 0.99. The inset shows the angular dependence of 
the SCE inelastic cross section for L = 8 [thin lines, Eq. ( I43b l and 
L ^ oo [thick lines, Eq. (EH, for different values of the incoming 
energy and interaction strength, as indicated in the legend. 



EIG. 5. (Color online) Inelastic scattering cross section for differ¬ 
ent incoming energies Em at fixed interaction strength U/J = 10 
and m = M, for a system of L = 8 sites at iz = 1, ver¬ 
sus scattering angle. Symbols correspond to exact numerical re¬ 
sults, while lines show the SCE analytical expression when the en¬ 
ergy gap En — Eq is evaluated to first order in J [Eq. (El, solid 
lines], and to second order [see Eq. ( lA4b . dashed lines]. We note 
that the maximum excitation energy in the system is max(i?„ — 
Eq) = 1.84i5r, and the number of contributing states (% of the 
total 6435 states) to the exact numerical cross section depends on 
the value oiEm/Er'. 2.73(100%), 0.6825(88.2%), 0.273(25.3%), 
0.1365(3.4%), 0.06825(0.4%). The analytical approximation (gill 
for L = 8 uses at the most 15 excited states \4‘q,x) to describe the 
cross section. 


independent value as L is increased. 

Here, we focus on the regime where the probe energy is 
high as compared to that part of the excitation spectrum of 
the system which contributes significantly to the cross sec¬ 
tion. Then, the corresponding excited states have comparable 
weights in the scattering signal [cf. the radicands in Eq. (fT^ l. 
Using the excitation gap to first order in J, as considered in the 
inelastic cross section ( |43]) . the high incoming energy condi¬ 
tion relevant for strong interaction reads i?in U{l — JEq]l:), 
which translates into 

y +2y/Mz/+l) + l. (71) 

In this regime the transferred momentum (1451) can be approx¬ 
imated by K ft! Kei, and, most interestingly, the inelastic cross 
section (1431) converges to a simple system-size independent 
expression as L ^ oo (see AppendixlBl): 


LT 

u ■ 

(72) 

The decay of the exact numerical cross section is compared 
to the latter analytical formula for different system sizes in 
Fig. |4l For strong interaction, Eq. (1721) describes remarkably 
well the numerical results, which converge very quickly to 
the L ^ oo cross section. We also note that for a typical 
choice of Ei„ ~ Ej., condition ( ItTI i translates into a regime for 
the interaction strength within which the L ^ oo expression 
should be valid; for jz = lit reads 1 U / J 148. Equation 
dT^ also reproduces correctly the overall angular dependence 
of the inelastic cross section even for not so high interaction, 
as can be seen in the inset of Fig. |4l 


1 da 


Nai dn 


inel 


= 8(jz + l)sin2 \ 


The existence of a system-size independent expression for 
the decay in the regime of strong interaction has a major con¬ 
sequence: Since Eq. (172] ! is able to describe the behavior of 
the inelastic cross section in the Mott insulating regime, we 
can conclude that, for incoming energies higher than the Mott 
gap, the cross section is non-zero inside the entire Mott phase 
in the thermodynamic limit, and does not vanish at the critical 
point. 

Up to second order in J, the analytical expression ( |4^ for 
the inelastic cross section is complete. It contains, however, 
higher incomplete orders in J, since the full functional de¬ 
pendence of the weighting factors — {En — Eo)/Ein, and 
correspondingly of the transferred momentum k, has been 
kept, where we simply inserted the expansion of the energy 
gap to leading order in J. These incomplete orders are only 
relevant for low incoming energy. An increasing deviation be¬ 
tween Eq. ( l43l l and the exact results should then be expected 
as Ein is reduced for a fixed interaction strength, as confirmed 
in Fig. |5] In this regime, an expansion of the energy gap to 
second order induces very little improvement, as also seen in 
Fig.El A higher order expansion of the matrix elements (|2^ 
would be necessary to describe the cross section more accu¬ 
rately for low incoming energy. 


C. Description for all 17/J: Bogoliubov result and SCE 

In Ref. we derived an analytical formula for the in¬ 
elastic cross section in the regime of weak interaction (small 
condensate depletion), making use of the Bogoliubov formal¬ 
ism. In combination with the SCE expression obtained in our 
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U/J 


FIG. 6. (Color online) Decay of the inelastic scattering cross section 
(logarithmic scale), as a function of the interaction strength U/J, 
for a system of L = 5 sites at filling factors v = 2 (upper panel) 
and V = b (lower panel). Symbols show exact numerical results, 
while lines correspond to analytical approximations: the Bogoliubov 
formula for small condensate depletion, Eq. (44) in Ref. 16911 (solid), 
and the SCE expression l l43t (dashed). Vertical dotted lines mark the 
position of the estimated critical points of the SF-MI phase transition 
tAppendix IaI. Relevant parameters are Ein = 2Er, 9 = 0.99 and 
m = M. 


A“"(dtr/dfl)i„ei [%] 

0 10 20 30 


5 
4 
V 3 
2 


1 

30 20 10 0 

U/J 

FIG. 7. (Color online) Minimum (dcr/dD).^^j of the relative dif¬ 
ferences (absolute value) between the exact numerical cross section 
and the analytical approximations: the SCE result [Eq. | |43H and the 
Bogoliubov expression [Eq. (44) in Ref. (H], as a function of the 
interaction strength U/J, for L = 5 and different integer bosonic 
densities. Eor every value of U/ J the relative difference is computed 
for the SCE expression [see Eq. I l70b 1 and for the Bogoliubov result, 
and then the minimum of both is taken. Symbols indicate the criti¬ 
cal value of the SF-MI phase transition, as predicted from third-order 
SCE (circles), and from DMRG calculations (triangles, for = 1, 2; 
data from Ref. l9^ ). Relevant parameters are E„, = 2Er, 9 — 0.99 
and m = M. 



present contribution, we provide an analytical description of 
the decay of the inelastic cross section which works remark¬ 
ably well almost in the entire range of interaction strengths, 
for any integer filling factor, as shown in Fig.|6l The shift of 
the regime of validity of Eq. (|4^ to larger U/J for higher v is 
compensated by the extended range of validity of the Bogoli¬ 
ubov approximation at high densities, Eq. (44) in Ref. 

The region where none of the analytical expressions seems to 
reproduce accurately the exact cross section can be correlated 
with the position of the SE-MI phase transition (in Eig.|6l the 
critical points estimated from SCE are indicated by vertical 
dotted lines). This can be clearly seen in Eig.|7] which shows 
A'^“((icr/dn)inei, the minimum of the relative differences be¬ 
tween the exact inelastic cross section and both analytical for¬ 
mulas, Eq. (44) in Ref. (Bogoliubov) and Eq. (|4^ (SCE), 
as a function of the interaction strength, for different filling 
factors. The approximate expressions deviate from the exact 
results in the vicinity of the transition, at which perturbative 
approaches can be expected to fail. Nevertheless, the overall 
analytical description performs remarkably well. 


V. CONCLUSIONS 

We have studied the inelastic cross section of a matter-wave 
scattered from a collection of interacting ultracold bosons in 
an optical lattice, focusing on the regime of large boson-boson 
interaction strength. We have employed a strong-coupling 
expansion (SCE) and a site-decoupling mean-field (ME) ap¬ 
proach to analytically describe the system eigenstates and en¬ 
ergies, from which analytical expressions for the cross section 


were obtained. The ME approach incorrectly predicts a linear 
decay of the cross section with increasing interactions, and its 
vanishing as the phase boundary from the SE to the MI regime 
is crossed. In contrast, SCE predicts a quadratic decay of the 
cross section, which vanishes only as U/J oo. Moreover, 
we derived an L —oo description of the cross section, which 
reveals that the inelastic scattering signal does not vanish at 
the SE-MI phase transition in the thermodynamic limit, but 
decays smoothly throughout the Mott phase. These findings 
are confirmed by exact numerical simulations. 

Together with the results previously obtained in the Bogoli¬ 
ubov regime ll^ . we provide a remarkable analytical descrip¬ 
tion of the decay of the scattering signal over almost the entire 
range of interaction strength, which seems to fail only in the 
vicinity of the critical point. The question remains open, as 
to whether or not the cross section bears a quantifiable finger¬ 
print of the phase transition. As shown in Eig.[8] the exactly 
computed inelastic cross section behaves in a clearly distinc¬ 
tive way depending on the phase of the system. Eor incoming 
energies larger than the Mott gap, the inelastic cross section 
will be non-vanishing in the Mott phase, as demonstrated in 
this work. Nevertheless, the transition in the thermodynamic 
limit could manifest itself as a non-analyticity at the critical 
point, i.e. a discontinuity of one of the higher derivatives of 
the inelastic cross section. This issue is not revealed within 
the perturbative approaches pursued in this work, and remains 
to be clarified. 
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FIG. 8. (Color online) Numerically calculated inelastic cross section 
for a system with L = 5 versus interaction strength and hosonic 
density n = N/L, where 3 ^ 25 (note the inverted abscissa 

axis). Black circles indicate the critical value of the SF-MI phase 
transition, as predicted from third-order SCE. Relevant parameters 
are Ei„ = 2Er, 9 = 0.99 and m = M. 


which after a very lengthy calculation can be found to be 
= -2Lv{v + 1) + ( 61/2 

— cos{qd) cos(2(^q — qd) [l + (L — 1) cos(2>!rd)] 

1j 

+ ^ [2u{u + 1) (6 cos{qd) - 1) + l] 

oJ-j 

+ sin(W) sin — 1)] |4n(n + l).5q,o 

-^^^^cos[<,(L-2)] 

2 (t/^ — l) - 

-— - - - cos \C,q{L - 2) + 2qd\ 

Lj 

+ cos [C,q{L - 2) + qd] 

X [(1 - (5q,o)(l + 2cos((?d)) - (L - 3)(5,,o] |, 

(A4) 

valid for L ^ 3, and where C,q is defined in Eq. ( |36] |. 

For the lowest energy fixed-density excitation of K^h, cor¬ 
responding to g = 0 (Cg = 0) and xd = -k jL, the SCE yields 
the energy gap [see Eq. for the ground state energy ex¬ 
pansion] 
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Ai =1 — 2 cos(7r/L)-\/l -F Av{v + 1) J 

+ J2|l -F 2v{v + 1)[3 - 2cos(27r/E)] 

+ Tpr sin2(7r/L)[5z/(z/ -F 1) -F 2]'l -F 0{J^), (A5) 


Appendix A: Strong-coupling expansion of \4>g,>€) states and 
estimation of the fixed-density SF-MI critical point 


The perturbative expansion of the eigenstates \(j)q,it) of 
ATbh which converge as J 0 to the one-PH states \q, x) 
[Eq. (|35] |1 is given in Eqs. (IMI l. The first order correction to 
the state, \4>^q}>c) is written in Eq. (i40l i. in terms of 


I \ i„\ 

X |n) > 


(Al) 


I ""®) =1111 


(g, V|E|n) {n\V\q,x) 

- Kl) 


(A2) 


which is system-size dependent. The critical point of the 
SF-MI fixed-density phase transition can be estimated from 
the vanishing of the energy gap in the thermodynamic limit, 
Aoo = 0. The expansion of Aoo is more easily obtained in 
a non-number-conserving approach, where the fixed-density 
excitation gap follows by adding the excitation energies from 
the ground state |^mi) to defect states containing one extra 
particle or hole. For these defect states, the SCE approach on 
a system with periodic boundary conditions provides a size- 
independent e^ansion for the gap valid in the thermodynamic 
limit iFt^ . From Aoo up to third order in J, the esti¬ 

mation of the critical point at integer filling factor v is 
given by the smallest positive solution of 

0 = 1- 2{2v + 1) J + [l + 2v{v + 1)] J2 + 2v{u^ + 2)P, 

(A6) 


The corrections appearing in the SCE for the energy [Eq. ( I39] l1 
up to second order in J are given in Eq. OUl, and 


= E 


|n)^Z> 


{n\V\q,x) 
fph - fn 


(A3) 


which evaluates to = (v/? — 2)/3 = 0.215 and = 
1/8 = 0.125. The estimates for these two critical points from 

DMRG calculations are = 0.305 and .iP = 0.180 iH. 
A comparison of the full phase diagrams ^ vs U obtained 
from third order SCE and DMRG calculations can be found 
in Ref. 1^. 
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Appendix B: Derivation of the L ^ oo SCE inelastic cross 
section 


Thus, the inelastic SCE cross section [Eq. (|43]) 1 for large in¬ 
coming energies can be approximated by 


The limit of large incoming energy, to first order in J, is 
given by the condition 

£;;„>£/(! +(Bl) 
which from the particle-hole dispersion (IJTI i translates into 

U 


^ - -f 2^Jl + Av{v + l). (B2) 


In this regime, the transferred momentum [Eq. (05])] becomes 
independent of q and k: k ~ Ke\, and the sum over x = 
TTj'/Ld, / = 1,..., L — 1, in Eq. (l43T l can be readily evalu¬ 
ated: 


1 da 

SCE 

■ J' 

Naj dn 

inel 

u 


8(z.+ l)|VE(Kel)| 
qd 


L2 


E 


sin 


|S(Kei-g)r- (B4) 


In the limit of a large lattice, we first note that 

-^|S(Kei-g)p —)> ‘^'^S{K^i-q-Q), 
L L-^-oo a ^' 

Q 


(B5) 


where Q = 27r j /d, j G Z, is a reciprocal lattice vector. More¬ 
over, the sum over the centre-of-mass quasimomentum q can 
be approximated by an integral: 


1 


E 


_^ d 

L—foo ‘2'K 


f'2TT jd 


dq. 


(B6) 


E l^[9>^(/)]r =4Esin2 . 


j'^1 


(B3) 


Using the two limits above in Eq. (IB4b leads to the E —> oo 
expression for the inelastic SCE cross section, Eq. ( |72] |. 
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